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Abstract 

We review the topological quantum computation scheme of Das Sarma et al. from the per- 
spective of the conformal field theory for the two-dimensional critical Ising model. This 
scheme originally used the monodromy properties of the non-Abelian excitations in the 
Pfaffian quantum Hall state to construct elementary qubits and execute logical NOT on 
them. We extend the scheme of Das Sarma et al. by exploiting the explicit braiding trans- 
formations for the Pfaffian wave functions containing 4 and 6 quasiholes to implement, for 
the first time in this context, the single-qubit Hadamard and phase gates and the two-qubit 
Controlled-NOT gate over Pfaffian qubits in a topologically protected way. In more detail, 
we explicitly construct the unitary representations of the braid groups S$a,, and and 
use the elementary braid matrices to implement one-, two- and three-qubit gates. We also 
propose to construct a topologically protected Toffoli gate, in terms of a braid-group based 
Controlled-Controlled-Z gate precursor. Finally we discuss some difficulties arising in the 
embedding of the Clifford gates and address several important questions about topological 
quantum computation in general. 

Key words: Topological quantum computation, Conformal field theory, Non-Abelian 
statistics 
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1 Introduction 



In contrast to the three dimensional world, where we could only find bosons and 
fermions, the statistics of localized objects in two dimensions turned out to be 
much richer, and includes fractional or anyonic statistics [1]. In the simplest case 
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of Abelian fractional statistics the counter-clockwise exchange of two anyons mul- 
tiplies the many-body state by the statistical phase exp(i6), where 6/n might be 
a fractional number. The most important distinction between anyons on one side 
and bosons or fermions on the other is that the clockwise and counter-clockwise 
exchanges are significantly different so that the many-body wave functions belong 
to the representations spaces not of the permutation group but of the braid group 
[1,2]. When the dimension of the braid-group representation is bigger than 1 the 
corresponding (quasi)particles are called plektons [3,4,5], or non-Abelian anyons, 
and the exchange of two such anyons results in a non-trivial statistical matrix acting 
over degenerate space of many -body states. As the matrices representing different 
exchanges do not commute in general, this kind of quasiparticle statistics is called 
non-Abelian. 

The most promising two-dimensional system, in which non-Abelian statistics may 
eventually be observed seems to be the fractional quantum Hall (FQH) state in 
the second Landau level with filling factor v = 5/2, which is now routinely ob- 
served in ultra high-mobility samples [6,7]. Convincing analytical [8,9] and nu- 
merical [10,11] evidence suggest that this state is most likely in the universality 
class of the Pfaffian FQH state constructed by Moore and Read [12] using correla- 
tion functions of certain operators in an appropriate 1 + 1 dimensional conformal 
field theory (CFT) [13] including the Ising model [12,14]. Experimental tests of 
fractional statistics appear to be much more difficult than those which confirmed 
the fractional electric charge [15] of the FQH quasiparticles. However, it turns out 
that consequences of the presence of non-Abelian excitations might in fact be easier 
to observe than the Abelian fractional statistics itself despite the structural compli- 
cations. In recent theoretical work interesting proposals for detection of the non- 
Abelian statistics of the quasiparticles in the v = 5/2 FQH state [16,17,18,19] and 
in the v = 12/5 FQH state [20,21], which is expected to be in the universality class 
of the k = 3 parafermion Hall state [22,23], have been made. 

It is quite remarkable that in addition to its fundamental significance, the non- 
Abelian quantum statistics might become practically important for quantum com- 
putation [24]. Although the ideas behind quantum information processing are sim- 
ply based on the well-established fundamental postulates of the quantum theory, 
its exponentially growing computational power could not have been exploited so 
far due to the unavoidable effects of quantum noise and decoherence as a result of 
interaction of the qubits with their environment. Even the quantum error-correcting 
algorithms [24], which allow to use operations containing certain level of errors, 
could not help creating a quantum computer with more than a few qubits. In this 
context recently appeared the idea of topological quantum computation (TQC) 
[25,26,27]. Because the interactions leading to decoherence are presumably local 
we can try to avoid them by encoding quantum information non-locally, using some 
global e.g., topological characteristics of the system. This topological protection of 
qubit operations means that quantum information is inaccessible to local interac- 
tions, because they cannot distinguish between the computational basis states and 
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hence cannot lead to decoherence and noise [25,27,28,29]. That is why topologi- 
cal gates are believed to be exact operations, which might potentially allow for the 
topologically protected quantum computation. 

The FQH liquid is a perfect candidate for TQC because it possesses a number 
of topological properties which are universal, i.e., robust against the variations of 
the interactions details, and could be successfully described by rational confor- 
mal field (RCFT) theories [30,31,32], which capture the universality classes of the 
FQH edge excitations (see also [33] and Refs. therein). The main idea is to use 
the braiding matrices [3,4,27] corresponding to the exchanges of FQH quasiparti- 
cles, acting over a degenerate set of topological many-body states, to implement 
arbitrary unitary transformations [25,26,27,34]. Because the single qubit space is 
two-dimensional, we need degenerate spaces of quasiparticle correlation functions 
with dimension at least 2. While the Abelian FQH quasiparticles have degenerate 
spaces on non-trivial manifolds such as torus [25,35], such constructions are not ap- 
propriate for experimental realization in planar systems, such as the FQH liquids. 
Alternatively, by Abelian anyons, one could realize in the plane only diagonal gates, 
such as Controlled-phase gates [36] but cannot implement non-diagonal gates such 
as the Hadamard gate. In contrast, the non-Abelian FQH quasiparticles by defini- 
tion have degenerate spaces even in planar geometry, which explains why the non- 
Abelian anyons might be more appropriate for TQC. The only residual source of 
noise and decoherence is due to thermally activated quasiparticle-quasihole pairs, 
which might execute unwanted braidings. Fortunately, these processes are expo- 
nentially suppressed at low temperature by the bulk energy gap, which leads to 
astronomical precision of quantum information processing. 

In a recent paper [28] Das Sarma et al. proposed to use the expected non-Abelian 
statistics of the quasiparticles in the Pfaffian FQH state to construct an elementary 
qubit and execute a logical NOT gate on it. The NOT-gate construction is very im- 
portant for quantum computation because it underlies both the construction of the 
single-qubit gates and of the Controlled-NOT (CNOT) gate [24] . However the NOT 
gate is certainly not sufficient for universal computation, though it is reversible, be- 
cause the universal classical gate, NAND/NOT, built from NOT is irreversible [37], 
while all quantum gates must be reversible. Therefore, in the TQC scheme of Das 
Sarma et al., we need to implement the CNOT gate (or reversible XOR), which 
plays a central role in the universal quantum computer [24,37]. 

The FQH state at v = 5/2 has one big advantage with respect to TQC — it is the 
most stable state, i.e., the one with the highest bulk energy gap, among all FQH 
states in which non-Abelian quasiparticle statistics is expected to be realized. Re- 
cent measurements of the energy gap [6,7] suggest that ratio noise/signal is ex- 
pected to be of the order of 10~ 30 or even lower, which is an unprecedented number 
in the quantum computation field and might potentially allow for the construction 
of a truly scalable quantum computation platform. On the other hand, one serious 
drawback is that the quasiparticles braiding matrices cannot be used alone for uni- 
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versal quantum computation because, as we shall see later, the representation of 
the braid group turns out to be finite. This has to be compared with the state at 
v = 12/5, whose braiding matrices are expected to be universal, but whose energy 
gap is an order of magnitude lower than the V = 5/2 one, which increases dramat- 
ically the noise/signal ratio. 

In this paper we extend the TQC scheme of Das Sarma et al. [28], which was 
originally based on monodromy transformations of the Pfaffian wave functions, in 
such a way to construct by braiding the single-qubit Hadamard and phase gates as 
well as the two-qubit Controlled-Z and Controlled-NOT gates. These constructions 
are naturally topologically protected. In addition we investigate some possibilities 
for topologically protected realization of the three-qubit Toffoli gate. 

Summary of results: we review in Sect. 2 the TQC scheme of Das Sarma et al. 
and introduce in Sect. 2.4 the 4-quasihole wave functions of Ref. [28] which we 
shall use to derive the elementary exchange matrices that will represent single- 
qubit gates. In Sect. 3 we consider in more detail the read-out transformation for 
the qubit of Ref. [28] by using the analytic properties of the 4-quasiholes states 
and prove, under certain conditions, the conjecture of Ref. [28] that the state of the 
qubit does not change after transferring one quasihole from antidot 1 to antidot 2, 
which is crucial for the construction, initialization and manipulation of the Pfaf- 
fian qubits. When these conditions are not satisfied, the TQC scheme of Ref. [28] 
is going to fail. We prove in Sect. 3.3 the orthogonality of the 4-quasiholes states 
forming the computational basis, which is very important for their quantum dis- 
tinguishability. Deriving explicitly in Sect. 6 the complete set of exchange ma- 
trices for the 4-quasiholes states, which has been partially done in Ref. [38] and 
completely reproduced in Ref. [39] using the underlying quantum group structure 
for the parafermion quantum Hall states, we construct in Sect. 8.3 all single-qubit 
gates, except for the n/S gate, entirely in terms of quasihole braiding. Then, in 
Sect. 8 we propose a natural two-qubit construction in terms of 6-quasiholes states 
and obtain explicitly the exchange matrices for these states. In Sects. 8.2 and 8.4 
we implement the Controlled-Z and Controlled-NOT gates entirely in terms of 6- 
quasiholes braidings. To the best of the author's knowledge this is the first explicit 
construction of these gates in the Pfaffian state, which is exact and topologically 
protected. We also construct in Sects. 8.6 and 8.7 the Bravyi-Kitaev two-qubit gate 
g3 and the non-demolition charge measurement gate, respectively, in terms of 6- 
quasiholes braidings. While the above gates are not sufficient for universal quan- 
tum computation they are known to form a Clifford group [40,35], which plays an 
extremely important role in error-correcting algorithms and, in particular could be 
efficiently used in such applications as quantum teleportation and super-dense cod- 
ing [24]. Moreover, if the Clifford group is supplemented by the so-called magic 
states or noisy ancillas that could already be used for universal quantum compu- 
tation [40]. In addition to the Clifford- group gates, instead of using the n/S gate, 
we propose in Sect. 9 to implement in a topologically protected way the Toffoli 
gate in terms of the CNOT and the Controlled-S gate, or by a braid-group based 
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Controlled-Controlled-Z gate precursor, which would actually form a universal set 
of topological protected gates realized with Pfaffian qubits. It appears that there 
is an additional topological entanglement in the three-qubit systems defined by 8 
Pfaffian quasiholes, which leads to complications in the embedding of the one- and 
two-qubit gates into systems with three or more qubits. This phenomenon seems to 
be common for all topological quantum computation schemes based on the braid 
matrices of the non-Abelian FQH anyons. 



2 The TQC scheme of Das Sarma et al. 



The main idea of Ref. [28,29] is to use the wave functions of the Pfaffian FQH state 
with 4 quasiholes to form an elementary qubit. Then quantum gates can be executed 
by braiding some of the quasiholes (i.e., by counter-clockwise exchanges of quasi- 
holes) leading to unitary transformations in the qubit space. When the positions of 
the quasiholes are fixed these wave functions form a 2-dimensional space which 
could be used as the single-qubit space. In general, the wave functions containing 
In quasiholes with fixed positions form a linear space with dimension 2" _1 [38]. 
The main reason for the exponential increase of the space dimension is the non- 
Abelian statistics of the quasiholes, i.e, when two quasiholes are fused together 
(taken to the same point in the coordinate plane) the result contains more than one 
quasiparticle due to the chiral Ising model fusion rule [12] 

ax o = I+i//. (1) 



The exponential degeneracy of the 2n quasiholes spaces can be alternatively under- 
stood by interpreting the Pfaffian state as a p-wave superconductor of composite 
fermions [8], where the non-Abelian quasiholes are represented by half-quantum 
vortices, and their non-Abelian statistics follows from the existence of Majorana 
zero-modes in the vortex cores obtained as solutions of the Bogolubov-de Gennes 
equations [8,9,41]. 

One way to keep the positions of the quasiholes fixed is to introduce antidots 
[42] (lithographically defined potential hills expelling the FQH fluid and creat- 
ing a "hole" or "island" inside it) in the FQH liquid and localize the quasiholes 
there as shown on Fig. 1. The positions of the quasiholes are denoted by rj a , where 
a = 1, ... ,4, and we assume that the quasiholes with coordinates r]\ and 772 form 
our qubit, while 773 and 774 (not shown explicitly on Fig. 1) are used to measure and 
manipulate the qubit's state [28,43]. 
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Fig. 1. (Color online). A v = 5/2 Hall bar with two antidots, on which the quasiholes (de- 
noted symbolically by small empty circles) with coordinates r/i and r/ 2 comprising the qubit 
are localized, and 4 front gates, M, N, P and Q creating tunneling constrictions. Black ar- 
rows depict the edge states, while the green and red arrows denote two alternative tunneling 
channels, a direct one and such enclosing the two antidots. 

2.1 State of the qubit and its initialization 

In order to have a real TQC scheme we need to specialize the computational basis 
{|0), |1)}. The definition of these states is closely related to the way we could 
possibly measure them and to the way we prepare the system in a certain state. 
That is why we shall start by saying how we can initialize the qubit and measure it. 

To initialize the qubit we put charge 1/2 on antidot 1 [28]. This can be done by 
adding one quantum h/e of magnetic flux in the vicinity of antidot 1 (e.g., by a 
solenoid piercing the antidot). A detailed analysis of the tunneling situation in the 
stable strong-coupling regime in a Pfaffian antidot with Aharonov-Bohm flux can 
be found in Ref. [18]. The FQH liquid containing the antidot responds by local- 
izing a charged excitation on the antidot border carrying one flux quantum. Be- 
cause of the fundamental FQH effect relation between the magnetic flux and elec- 
tric chargeQ 

2ei = v4>, v = - 

the charge of such an excitation is 1/2. There are only two allowed charge 1/2 
excitations of the Pfaffian FQH state which could be localized on the antidot: in the 



1 as usual we consider only the fractional part of the filling factor v = 2 + 1/2 correspond- 
ing to the top-most Landau level 
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notation of Refs. [33,18] these are represented by the field operators 



Hz). 



and ty/(z):e'^ 



Hz). 



<£=!, <2el=lA 



(2) 



where is the Ising model Majorana fermion and the vertex exponent of the 
normalized boson field 0(z) represents the u(l) part of the excitation. Now we can 
envision the following computational basis 



or, in other words, the state of the qubit is |0) if the Majorana fermion is not oc- 
cupied and |1) if it is occupied. Note that the states corresponding to the fields in 
Eq. (2) cannot form coherent superpositions because they belong to different super- 
selection sectors inside the Neveu-Schwarz sector (when the flux threading antidot 
1 is one quantum in the experimental setup), which correspond to their different 
fermion parity. 

In Table 1 we give for convenience the list of the 6 topologically inequivalent sec- 
tors (quasiparticles) for the Pfaffian FQH state. Note that the topological sectors in 
which the chiral fermion parity is well-defined contain quasiparticles with both pos- 
itive and negative parities because they can be obtained from each other by adding 
an electron. The quasiparticle spectrum is obtained from the chiral partition func- 
tions of the corresponding topological sectors, see Refs. [33,18] for more details. 



2.2 Measurement of the qubit state 

It was one of the bright ideas of Ref. [28] to use the electronic Mach-Zehnder inter- 
ferometry [44,45,19] to determine the state of the qubit. More precisely, let us try 
to measure the diagonal conductance, a xx , which is proportional to the probability 
for a charged particle to enter the lower edge of the Hall bar in Fig. 1 and to exit 
in out of the upper edge. Because of the constrictions created by the front gates M 
and N, P and Q, there are two alternative channels for a charged quasiparticle en- 
tering from the lower edge to exit from the upper one: one is to tunnel between the 
front gates M and N and the other is to tunnel between P and Q. Therefore, to lead- 
ing order in the tunneling amplitudes tMN an d tpQ, which are assumed to be very 
small, the diagonal conductance would be proportional to the interference of the 
two amplitudes [28]. Consider, for example the interference of the two tunneling 
processes if the charge 1/2 on antidot 1 does not contain Majorana fermion. Then 
the amplitude for tunneling between P and Q must be multiplied by the Aharonov- 



|0> « 
|1> < 



if the charge 1/2 state is 
if the charge 1/2 state is 




(3) 
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Table 1 

Topologically inequivalent quasiparticles in the Pfaffian FQH state and their quantum num- 
bers: electric charge Q, neutral-sector chiral fermion parity Jf, conformal dimension A and 
quantum statistics 6 j% = 2A mod 2 



Particles in the same 


Fields 


Charge 


Parity 


CFT dim. 


Quant, stat. 


Topological Sector 




Q 


Yf 


A 


e/n 


vacuum 


I 





+ 








hole/electron 




±1 


— 


3/2 


1 


quasihole (vortex) 


e ! '^ a 


1/4 


undefined 


1/8 


1/4 


quasiparticle 


e"^ a 


-1/4 


undefined 


1/8 


1/4 


+ 1 flux 


e ^ 


1/2 


+ 


1/4 


1/2 


— 1 flux x Majorana 


e -Ji iff 


-1/2 




3/4 


-1/2 


-1 flux 


e 


-1/2 


+ 


1/4 


1/2 


+ 1 flux x Majorana 


e ^ \jf 


1/2 




3/4 


-1/2 


+2 flux (K"-boson) 




1 


+ 


1 





Majorana fermion 


¥ 







1/2 


1 



Bohm phase for the non-Abelian quasihole which when tunneling between P and 
Q actually encircles the charge 1/2 state on antidot 1, i.e., 



e'7 



(4) 



because a quasiparticle with charge Q e \ = 1/4 encircles magnetic flux 4> = 1. 
Therefore, the diagonal conductance, which is proportional to the modulus-square 
of the amplitude, reads 



Oxx \tMN + itpQ\ 2 . 



(5) 



If instead, the state of the qubit is |1), i.e., there is a Majorana fermion on anti- 
dot 1, then in addition to the Aharonov-Bohm phase (4) there would be a minus 
sign coming from the fact that the Ising model o filed is transported around the 
Majorana fermion. To see, in the CFT language, why this minus sign appears, we 
consider the operator-product-expansion 



o(0) 



Tz 



so 



that (7(z)y(0) — > — ct(z)v(0) when z 



■ e 2Ki z, 
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i.e., when o is transported adiabatically around i//. Thus, the diagonal conductance 
measurement in the state 1 1 ) gives 

Oxx x \tMN-itpQ\ ■ (6) 

Note that the two different interference patterns, Eqs. (5) and (6) of the diagonal 
conductance are very well distinguished experimentally due to the high visibility 
of the Mach-Zehnder interferometry [44] . 

2.3 Splitting the 1/2 charge: finalizing the qubit 

Despite that we can efficiently measure the two states in the computational ba- 
sis this is still not sufficient for quantum computation. The reason is that we need 
to form coherent superpositions of the states |0) and |1) which is not allowed for 
the charge 1/2 states (3) due to the fermion parity superselection rule. In order 
to circumvent this difficulty Das Sarma et al. have made another interesting pro- 
posal [28]: to split the charge 1/2 state into 1/4 x 1/4 state by transferring one 
charge 1/4 from antidot 1 to antidot 2. This is indeed possible, if one applies volt- 
age between the two antidots, because the most relevant quasiparticle for tunneling 
through the bulk of the Pfaffian FQH liquid carries charge 1/4. This quasiparticle is 
non-Abelian and contains a o field from the Ising model because there is no other 
charge 1/4 quasiparticle and the non-Abelian one has the minimal CFT dimension 
[12,38,33]. Now the state on antidots 1 and 2 is equivalent to o{r\\)c(r\2) which 
together with the quasiholes at 773 and 774 correspond to a 4-quasihole wave func- 
tion, which belongs to a two-dimensional space as discussed in the beginning of 
Sect. 2. 

While the charge 1/2 states (3) cannot form coherent superpositions, because they 
belong to different superselection sectors, for the 4 quasiholes configurations we 
can consider the superposition of the states obtained by splitting of the vacuum in 
two different ways, namely 

I^IxI^(axa)x((jxff) 

I^fXl(/->(ffXff)x((JX(T), (7) 

because the two 4-quasihole states now belong to the same superselection sectors. 

One natural question arises in this charge splitting procedure: does the state of the 
qubit remain the same in the process of transferring one charge 1/ '4 from antidot 1 
to antidot 2 or it changes? Or, even, more philosophically: does the system with 4 
O fields has something in common with the computational basis of our charge 1/2 
state plus 2 additional o fields, or it is completely different? 
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The original idea of Ref. [28] was that the state of the qubit does not change during 
the splitting procedure because the pairs of quasiholes in the Pfaffian state share 
a pair of Majorana fermions zero modes, (l//b, ty/b), whose combined (left- plus 
right- moving) fermion parity is supposed to be conserved. That's why the read- 
out should give the same results as before splitting unless some Majorana fermions 
could tunnel from the edge or from another antidot, however, these processes are 
exponentially suppressed because the quasiholes are supposed to be well separated 
and far from the edges. 

In Sect. 3 we shall explicitly derive the readout results for the 4-quasiholes wave 
functions and will find conditions under which this is in agreement with Ref. [28]. 

2.4 Four-quasiholes wave functions 

As we shall see in Sect. 3, the topological phase, which determines the conductance 
interference pattern, can be obtained from the monodromy matrices appearing in 
the adiabatic transport of some quasiholes around others, so that we need the wave 
functions explicitly. The wave function for even number N of holes (or electrons) at 
positions z\ , . . . , Zn containing 4 quasiholes at positions r/i , . . . , 774, can be realized 
as a correlation function 

N 

^4qh(T?l, T? 2 , 13, Ui}) = {Vqh(f7l)Vqh(m) VqhMVqhfo*) Vholefe)) 

i=l 

of the field operators corresponding to creation of holes and quasiholes 

VWz) = V^):e iV5 ^: and yftft) = aft) : e^ 0(l?) :, (8) 

respectively, where o(r]) is the chiral spin field in the Ising model of dimension 
1/16 and y{z) is the right-moving Majorana fermion in the chiral Ising model. It 
can be expressed as 



X i , 4qh(rihm, r l3,n4;Z h ...,ZN) = ( C7(T7i)(7(T]2)<7(T]3)C7(T74)n^;) ) x 

I N 

n ni nvfc-n) (zi - vi) 5 - t? 3 ) 5 - 774) n ^ 

l<a<b<4 i=l l<i<j<N 

where the average sign (• • ■) now represents the vacuum expectation value in the 
chiral Ising model and the last three product factors in Eq. (9) come from the w(l) 
components of the holes : exp(?'v^20(z)) : an d of the quasiholes : exp(/0(z) /2\^2) :. 
We used here the notation r\ a i ? = r\ a — r\i, for a ^ b and zij = z\ — Zj for i ^ j. 
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One important detail is that the chiral field o(rj) does not have a definite fermion 
parity because of the non-Abelian fusion rule (1) which mixes states with different 
fermion parities. Nevertheless it would be convenient to introduce two chiral fields 
<7± of dimension A± = 1/16 and opposite parities [33] 

7fo±(t7)7f = ±o±(t7), (10) 



in terms of which the non-Abelian fusion rule (1) necessarily splits into two Abelian 
channels 

<y± x o± = I, o± x a T = yr. (11) 



The a filed entering Eq. (9) is then identified with 



In order to obtain the 4-qh wave function (9) we carefully repeat the arguments 
of Ref. [38] in which the chiral many-body wave functions of the Pfaffian state are 
obtained by bosonization techniques from the c = 1 complex Ising model. The final 
result is 

^4qh(m, 12, 13, l4-,Z h . ..,Z N )= ^Sh + ^l ( 13 ) 



where we have used the notation of Refs. [28,38] = |0), *¥%L = 1 1) would be 
our computational basis in the 4-qh wave function's space) 



with x being a CFT invariant crossratio [13] 

X= ^^ and (lab = la- lb), (15) 
^13^24 



^ )( «o = Pf ^ — J x 

x n fe-^-) 2 5 

l<i<j<N 
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where {a,b,c,d} is a permutation of {1,2,3,4} satisfying a < b and c < d. The 
Pfaffian of an anti- symmetric matrix Mu of even dimension N is defined as 



1 



N/2 



Pf(My) 



2 N / 2 ((N/2)\) 2 



£ sign(cr)f[M CT( 2 fc _i) CT (2fe). 



It is worth stressing that the space of the 4-quasihole wave functions for fixed posi- 
tions of the quasiholes is two-dimensional. The second independent wave function 
could be obtained from the first one, Eq. (13), by transporting 773 around 774 (i.e., 



T734 e 27r; T] 3 4) which transforms ^ U Q j h -> ^ J h and -> -*P 4 ' j h so that 



It appears, however, that in some situations, such as the qubit initialization of 
Ref. [28], the 4-qh wave function may be driven directly in the state W^tL or 'Pjj!!. 



3 Read-out in the TQC scheme of Das Sarma et al: the measurement of the 
4- quasiholes qubit state 

As we have seen in Sect. (2.2) the read-out of the qubit is performed by interference 
measurement of the diagonal conductance, for which we need the monodromies of 
the corresponding wave functions. In this section we shall compute the explicit 
monodromies of the 4-qh wave functions (14) as well as for (13) and (16). Then, in 
Sect. 4, we shall compute the corresponding monodromies for the two charge 1/2 
states, which can be obtained by fusing quasiholes with coordinates r\\ and 772. 

Assuming that our qubit is formed by the quasiholes with coordinates r\\ and 772 
we can interpret the quasihole with coordinate 773 tunneling either through M and 
N or through P and Q as generating the interference pattern in the longitudinal 
conductance, like in the non-Abelian Mach-Zehnder interferometer [45]. In more 
detail, if the third quasihole tunnels through M and N, as shown on Fig. 2, this 
could be interpreted as a clockwise braid of this quasihole with the qubit, while if it 
tunnels through P and Q this gives rise to a counter-clockwise braid (denoted below 
as R; cf. Ref. [45]) so that the quantum amplitude for these two processes is (we 
have absorbed all dynamical phase factors for the two paths in the corresponding 
amplitudes tMN an d tpQ) 



*4qh(m,f]2, *?3, r]4-,Z h . ..,Z N )= rf h ~^ 



•(i) 

4qh" 



(16) 



\A)=t MN R 1(^1,772), r/ 3 , 774) +^^1(771, 772), 773, 774) 



12 



► ' 

Fig. 2. Conductance interference term in the read-out procedure expressed by the braids R 
(corresponding to tunneling through P and Q), R 1 (corresponding to tunneling through M 
and N) and the monodromy R 2 . 

and the longitudinal conductance would be 

o xx oc( A \A) = \t MN \ 2 + \t PQ \ 2 + 2Rc (t* MN t PQ (y\R 2 \y)) , 

where = | (rji, 772) , 773, 774) is the 4-qh state in which the qubit is formed by the 
quasiholes with positions 77 \ and 772, and we have used the unitarity of the braid 
operator R and R 2 is the corresponding monodromy operator. The operator R 2 is 
actually the operator which takes the quasihole with coordinate 773 around those 
with r/i and 772 in counter-clockwise direction. Note that if we fuse the quasiholes 
forming the qubit, i.e., 771 — > 772, and the state is I then the matrix element is +i, 
while if they fuse to i// the matrix element is — i, which reproduce the interference 
results in Sect. 2.2. It is worth stressing that there is a remarkable relation between 
the expectation value of the monodromy operator,/? 2 , corresponding to the adiabatic 
transport of particle with label a around particle with label b, and the modular S 
matrix [13] associated with any rational CFT [20]: 

(ni\R 2 \\fi\ = SabS °° 

SoaSob 

where the label corresponds to the vacuum. Working with modular S matrices 
is very convenient because they are explicitly known for almost all rational CFT 
related to the FQH effect. 

Thus we see that the read-out or measurement procedure for the 4-qh wave func- 
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tions consists in taking the quasihole with coordinate 773 around the qubit [28,43], 
i.e., around the two quasiholes with coordinates 771 and 772. E] 

It is worth stressing that the 4-qh wave function basis (14) is very convenient from 
the point of view of the adiabatic transport because the functions ^(ab)( c d) nave no 
Berry phases [38] so that the transport effects could be explicitly determined from 
the monodromies of the multivalued function 

f{z) = y/l±vk, zGC, 



entering the denominator in Eq. (14). Notice, however, that "taking 773 around 771 
and 772" in the functions (19) for z = x with x defined in Eq. (15), which are simply 
proportional to the 4-point functions of the chiral spin field in the Ising model, is 
equivalent to the adiabatic transport of 773 around 774 alone. This is because the 
monodromy for the transport of 773 around 771, 772 and 773 is trivial, as shown on 
Fig. 3, since the contour on the left can be contracted to a point (at infinity) without 




Fig. 3. The trivial monodromy of the 4-quasihole wave function, when one quasihole is 
transported around all others, allows to compute the read-out as taking 773 around 774. 

passing through any other singularity. Therefore it would be simpler to compute the 
read-out transformation of the functions (19) for z = x by 

R 2 34 - 7734^ 77^4 = f lime 27r! ' f 7734, tE[0,l). (17) 



The transformation (17) obviously preserves the absolute value of the crossratio 
x = 77127734 = \-x <£> x=l-x, (18) 

f?13T?24 



where x is defined in Eq. (15), because all other 77^ = r\ a b, i.e., the contour for 
transportation of x, corresponding to Eq. (17), is a circle with center at x = and 

2 Alternatively, as a matter of choice, one could use the quasihole with the coordinate 774 to 
encircle the quasiholes localized on antidots 1 and 2, which would lead to the same results, 
however, we will stick here to the notation of Ref. [28]). 
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radius |x|. Therefore the result crucially depends on whether |x| is bigger or smaller 
than 1 . Thus we need to consider the behavior of the functions 

f ± (x) = \Jl±VT^Z (19) 

under the transformation (17). Yet, it is more convenient to first analyze Eq. (19) as 
a function of z and then just outline what happens when we change z = l—x. 

The multivalued function (19) has two separate branching points: one is z = 
which, if encircled by the continuation contour, would change the sign of the inner 
root, as can be seen from Eq. (A. 2) in Appendix A, i.e., 

yfl±Vz f o r z^e 2ni z and | z | < 1. 

However, the second function, corresponding to the minus sign under the square 
root in Eq. (19), has one more branching point at z = 1 which, if encircled in the 
process of the analytic continuation changes the sign of the outer root, i.e., 

y/l±Vz -> ±\Zl±Vz, for (z-l)^e 27r/ (z-l) and|z-l| < 1. 

This can be easily seen from the Laurent-mode expansion for \z — 1 1 < 1 , by looking 
at Eq. (A.4) in Appendix A. 

Now consider what happens when we change z= l—x. When x — > e 2mt x, and t 
goes from to 1, this transports x along a circle of radius |x| and center at x = 0. 
Then — x is transported along the same circle, though with a phase shift of %. Thus 
z= 1 — Jc is transported along a circle of the same radius however translated to 
be centered at z = 1 as shown on Figures 4 and 5. 

Therefore we need to consider both cases, |x| < 1 and |x| > 1, separately in more 
detail. Of course, one may consider other contours which are homotopic to that 
of Eq. (17) but do not preserve the absolute value of x. In that case the results 
would be the same as those obtained from Eq. (17), only the conditions on |jc| 
must be replaced by the homotopic condition whether these contours encircle both 
branching points z = and z = 1 or only one of them. 



3.1 The read-out for \x\ < 1: take T73 around 7)4 

According to our analysis, when |x| < 1, the read-out transformation (17) corre- 
sponds to transporting the crossratio (18) along a contour which encircles only the 
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branching point at z = 1 , 

Cf = {z=l-Q lKit x\0<t< 1, |x|<l}, (20) 

as shown in Fig. 4. We have introduced two branch-cuts for the function (19) like 













Ixkl 


branch-cut 










\ branch-cut 


(—,0) 


z=0 






z=l 


1 OL+oo) 



Fig. 4. The contour used for the analytic continuation of Eq. (19) when |x| < 1. As 
it encloses only the branching point at z = 1, going along this contour only changes the 
outer-root sign of 

in Fig. 4. Because the 4-quasihole wave functions v f(i3)(24) and v I / (i4)(23) have no 
branch-cuts in T]^, they acquire no phase under the transformation (17), and the 
only phases come from the prefactor in Eq. (14) containing rj a b and x. 

Since |x| < 1 and the transformation (17) transports x along the contour (20), which 
is shown on Fig. 4, this transformation only changes the outer-roots signs, i.e., 

p2 f Vi + VT^ \ _ (I 0\ f v /i + ^Y^ \ 



m(0) vp(l)' 
T 4qh' T 4qh 



we 



In order to find the action of the transformation (17) in the basis 

only have to add the phase coming from T]^ 4 in Eq. (14) which is & l7t / 2 . Notice that 
we use "taking 773 around 774" as the readout prescription only for the functions 
(19) because these are the functions for which the total monodromy is trivial. For 
all other functions, including fractional powers of T\ a b, we still use as the read-out 
"taking T73 around r\\ and 772". Thus we obtain the read-out transformation as 



(0) 



TJ i\x\<\) I ^4qh 
read— out | > T/ (1) 



4qh, 




for 



7?12T?34 



r?13?724 



< 1. 



(21) 
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Notice that this "±i", on the diagonal of the matrix in Eq. (21), is exactly the topo- 
logical phase that appears for the quasiparticles traveling along the contour passing 
through the front contacts P and Q, see Fig. 1, in the diagonal conductance mea- 
surement so that 

<yS <l « \tMN±itpQ\ 2 , with " + " for |0) and "-"for |1). (22) 
3.2 The case \x\ > 1 

Alternatively, if |x| > 1, then the read-out transformation (17) leads to transporting 
x along the following contour 

= { z =l-e 27tit x\0<t< 1, |x|>l}, (23) 

which encircles both branching points at z = 1 and z = as shown on Fig. 5. 
Therefore, the read-out transformations (17) would change the sign of the outer 




Fig. 5. The contour C~ used for the analytic continuation of Eq. (19) when |jc| > 1. Because 
it encloses both branching point at z = 1 and z = going along this contour changes the 
outer-root sign whenever the inner one is "— " and flips the inner-root sign of f±(z)- 

root whenever the inner sign is "— " as well as would flip the inner-roots signs, i.e., 
D 2 fy/ l + y/l-x \ _ ( l\ / Vl + Vl-x ^ 

R ^{V^7T=x)-{-i o){^7T=x)- 
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Adding again the Abelian phase e l7t / 2 , coming from r\\!^ in Eq. (14) we obtain the 
read-out transformation as 





T?12T?34 



T?13T?24 



> 1. (24) 



Now the diagonal conductance measurement might be different from Eq. (22), as it 

crucially depends on the overlap between the states (14). If the two states are 
orthogonal, which is a fundamental requirement in quantum theory, then the mon- 

odromy average (^4qh I -^34 1 ^4qh ) ' determining the interference pattern, would 
vanish and the diagonal conductance could not distinguish between the two states 
in the computational basis. Thus we see that the Mach-Zehnder interference mea- 
surement could only work for |rj 12 T734 1 < 1 7713 7724 1- 

Using the CFT invariance of x, one can prove that the absolute value of the cross- 
ratio x depends on the absolute values of the quasiholes positions, i.e., 

|x| < 1 <^> 1 773 1 < |t7 2 |, while |jt| > 1 <^ 1 773 1 > |T7 2 |. 

Remark 1 The main result in this Section has very important implications for the 
experimental realization of the topological quantum computer in FQH systems at 
v = 5/2. The read-out procedure crucially depends on the absolute value of the 
crossratio (18). When \x\ < 1 the TQC scheme of Ref. [28] is going to work as 
originally proposed, while if \x\ > 1 this scheme would fail. As the absolute value 
of x depends on the quasiholes positions, this certainly gives some hints on how 
the antidots should be arranged in order for the read-out procedure to work as in 
Ref. [28] and, as we shall see later, this is also related to whether the NOT gate 
could be executed in the Pfaffian qubit or not. 



3.3 Orthogonality of the 4-quasiholes wave function 



One of the fundamental requirement in quantum theory is that the states |0), |1) 
are not only linearly independent but in fact orthogonal. Otherwise there exists no 
measurement that can reliably distinguish these states [24]. Here we will demon- 
strate that the 4-quasiholes states (14) are indeed orthogonal, at least for |x| < 1. 
This seems natural because, as we already know the Mach-Zehnder interferometer 
can distinguish the states (14), hence these states should be orthogonal. However, 
it is still worth verifying this directly if possible. 

The point is that for |x| < 1 we have, according to Eq. (21), U^^fl = and 
u¥ 4 l = , hence, U^[% = -i¥ 4 % and U*P® = where U = 
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^read-out- In order to find the overlap ^4^h»^4qh) we ^ rst P^ u § i ns ide ^ me oper- 
ator U and use the above relations, i.e., on one side we have 

while, on the other side, using the definition of the Hermitean conjugation and the 
anti-linearity of the inner product with respect to its first argument, it is equal to 

(^IXL) = («) = ■ <*> 

Obviously the right-hand sides of Eqs. (25) and (26) could only be equal if 

^ T 4qh' T 4qh J ~ V ' 



4 Measurement of the charge 1/2 state's wave function 



As we have seen in Sect. 2.2 the state of the qubit before splitting 1/2 —> 1/4 x 1/4 
is |0) if the Majorana fermion is absent or 1 1) if it is present in the charge- 1 /2 state. 
Using the fusion rules for two quasiholes (8) 

Mm) Vqh(r? 2 ) ^ (l + J^Yim)) e'^ Hm) (2V) 



as well as the general fusion rule of vertex exponents 

e j\<Kn a ) e *'4<Kf?i>) hJ{X a +h)<t>(nb) 



we can deduce the form of the charge 1/2 wave functions from the 4-quasihole 

wave function (9) by taking r)\ — > 772. Let us denote by and *Fi the result of 

.j>_ .j>_ 

fusing the two quasiholes \j/ qh to e' ^ and \j/e ^ for T]i — > r\2 as in Eq. (27), respec- 
tively, i.e., 



W = ( ofoWrriflYfaM 



1 mo 1 -iL \ ^ <T,4) A • /o*/ 



!=1 / \ (=1 
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Next we use the well-known formula 

ia=l / a<b 



to show that 



A iJiAhA . 1/4 . 1/4 . 1/8 /independent of' 



e*H* (,b) e* (,h) ne^ \ = ^{W^ 

i=l I \ T]ab 



where the second expectation value on the right-hand side is independent of r\ a b and 
therefore gives no contribution to monodromies. Similarly, noting in passing that 
the branch-cut structure of the wave functions with 2 or 3 quasiholes is determined 
by the corresponding 2-pt and 3-pt quasihole functions, for even N we find 



_i/8 / single — valued \ 

'=1 / \ mriab I 



t—w \ —1/8 / M11 & lc — VcllUCU \ 

CT(TJ 3 )a(TJ 4 )n^«))=T?34 ( > and 



A \ r\l ! * /single- valued \ 

Kl2)ff(l3)ff(l4)l]^i) ) = 1/2 1/2 ' 



I 



' =1 / ^23 ^24 \ m Vab 

Thus we finally obtain, for even number /V of electrons, 



i / single — valued \ 

X I , 0=(T?23T?24)M ) 
\ ln lab I 

y = V^hA I single - valued \ 
(l?23»?24)M inT U / 

We point out that the factors containing Tl a / 7 , hence the monodromies of the wave 
functions (29), are independent of whether the absolute value of x is bigger or 
smaller than 1 because considering 7712 — >■ actually means |jc| < 1. 

Now the read-out procedure, Eq. (17), which after the fusion r/i — > T72 is reduced 
to 7723 — ► Q 2m T}23, simply gives for the channel passing through P and Q in Fig. 1 

¥0 -> e'^o, while -> e^^i 
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so that the diagonal conductance is in agreement with Eqs. (5) and (6) 

Oxx < x \tMN + i tpQ 1 2 for , while o xx <x\t MN - i t PQ \ 2 for ¥ i . (30) 



5 The NOT gate of Das Sarma et al. 

The original idea of Ref. [28,29] behind the construction of the logical NOT gate in 
the setup shown on Fig. 1, is that when the quasiholes with coordinates 7]i and 772, 
localized on antidots 1 and 2, form the qubit [43] another quasihole with coordinate 
773 could tunnel between the front gates A and B (through an additional antidot 
located between A and B in order to guarantee a single quasihole tunneling) in 
such a way to flip the state of the qubit. The point is that this third quasihole actually 
encircles the quasihole localized on antidot 1 but not the second quasihole localized 
on antidot 2. More generally, while the read-out of the qubit is performed by taking 
r/3 around 771 and 772, the NOT gate could be executed by taking T73 around 771 or 
772 but not both of them. 

To understand this in more detail let us consider the simplest situation when the 
quasiholes are ordered \ r]i \ > \r\2\ > | T73 j > | T74 1 and the quasihole at 773 traverses a 
closed loop around the quasihole at 772, i.e., 

NOT: 7723 lim e 2Kit 7723, where fe[0,l) (31) 

with all other quasiholes coordinates remaining unchanged. Provided that the three 
quasiholes at 771, 772 and 773 are kept well- separated, after the tunneling of the third 
quasihole the 4-qh wave functions (14) are adiabatically transformed into new wave 
functions which we shall find now. First the transformation (31) transforms the 
crossratio (15) according to x — > e x. Next we note that the Pfaffian wave func- 
tions ^(afytcd) are single-valued, so that they remain unchanged under the NOT 
transformation (31). Finally, the square root of x changes sign under the NOT trans- 
formation (31), i.e., 

s/x -»• Ve 2 ™x = V*. 

The transformation of a/1 ± yfx is more subtle and depends on the absolute value 
of the crossratio x. When |jc| < 1 we have yl± y/x — > a/1 =F y/x un der the transfor- 
mation (31), while if |x| > 1 \J\ ± yfx — > ± a/ 1 =F y/x, because the transformation 
contour now encircles both branching points. Thus we find that the two 4-qh wave 
functions (14) transform under the NOT transformation (31) as follows 
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not (< , r vi ;^_ 



^"-(jIM-i o)C|) f ° r w> l <33) 



Thus we conclude that the transformation (31) indeed maps <-+ ^qh- 



6 Exchange matrices and monodromy group representation of the Ising 4- 
point functions 



The four-point correlation functions of the chiral spin field (12) in the Ising model 
can be shown to be [46] 



F(VhV2, t?3, r? 4 ) = <<t(77i)ct(t72)<t(77 3 )o-(t74)) = F+ + where 

are the chiral conformal blocks, which could be expressed in terms of the fields o± 
with definite fermion parity as F± = (a + (rji)a±(r/2)cr + (773)o"±(774)). 

The counter-clockwise exchange of the quasiparticles with coordinates T]i and T]2 
can be preformed by analytic continuation along the contour defined by 

^mpi+^mzik, r& = mp>-&mzi!i, o<r<i,<35) 



as shown in Fig. 6. Executing this transformation and taking the limit t — > 1 we 
obtain for the 4-pt functions (note that r\[ 2 = e i7r r]i2) 



^12^14^23^34/ V 13^24 / V V ^ 14^23 



Therefore (fixing the signs of the square roots by the requirement that the i?-matrices 
must satisfy the Yang-Baxter equations (37)) we get a diagonal matrix 



F V i \F- ' 



Precisely in the same way we can compute the exchange matrix R34 which takes 
precisely the same form in this basis. In order to compute the exchange matrix R23 



22 




Fig. 6. Counter-clockwise exchange R12 of coordinates T]i and r/2 performed by taking the 
limit t = 0^ 1_ inEq. (35). 

we apply the coordinate transformation (35), after renaming (771,772) — > ( r ?2, T 73) 5 
and use the identity 



1 + \fx + iX y 1 — \fx = vlXiy y/l — x — iX \/x 



for A = ±1 . We obtain 



Ri3F± 



T713T724 



y/2 \ T?1277l47?23 7734 



v/2 



e ±; 4 (f + q= ?F 



Thus, we can summarize our results for the elementary exchange matrices i? flia +i 
in the basis {F + ,F_} 



^23 



e J s / l -i 



V2\-i 1 



(36) 



All other exchanges can be obtained from the elementary ones, e.g., 

it 

e 8 
71 



_i e ! s / 1 1 



etc. 



(0,1) 



f5. i Exchange matrices for the 4-quasiholes wave functions 



The 4-quasiholes wave functions (14) are built up from the 4-point functions of 
the Ising model and the functions ^{ab){cd) which are single-valued in the positions 
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of the quasiholes. Now a natural question arises: is the braid-group representation 
for the 4-point functions (19) extended to a braid-group representation over the 
4-quasihole functions? The answer is yes. The point is that when a quasihole tra- 
verses a closed loop, the 4-quasihole functions should acquire an additional phase 
proportional to the number of electrons inside the loop (or, to the area of the loop). 
This might lead to a projective representation of the braid group as looks to be 
the case in the the p-wave composite-fermion superconductor approach, where the 
quasihole is identified with a half-quantum vortex [9]. However, as can be seen 
in our approach directly from the 4-quasihole functions (14), each electron inside 
the loop contributes 2n to this phase because the quasiholes are by definition lo- 
cal with the electrons^. Thus the phase is insensitive to the number of electrons 
inside the loop, it only counts the number of quasiholes. Nevertheless it turns out 
that braid-group (or mapping class group) [2] representations in terms of CFT cor- 
relation functions are generically projective. The point is that the coordinates of the 
many-body wave functions, which in the CFT approach [33] are chiral correlation 
functions defined on the unit circle, could be naturally extended by analytic contin- 
uation to the vicinity of the unit disk. Then by CFT transformation these functions 
could be extended to the entire compactified complex plane, which is isomorphic 
to the two-dimensional sphere. Now, besides the Artin relations [2] 

B i B J =B j B i , for \i-j\>2 

B i B i+1 B i = B i+1 B i B i+h where B,=%iG^, (37) 

for the generators B; (i = 1, . . . ,n — 1) of 3$ n , the representation of the braid group 
on the sphere should satisfy one more relation [2] 

B\B2 - B n -2B 2 n \B n 2- • B2B1 = I. 

As can be seen by direct computation the above relation is satisfied by the elemen- 
tary exchange matrices of the chiral CFT correlators only up to phase so that the 
braid-group representation is projective. 

In the rest of this subsection we will show that from the exchange matrices for the 
functions F± we could obtain the corresponding matrices in the basis {^4^, ^4^} 
in the form [47] 

M 4 2=4 4 4=(o f). *23=^(-i l)' (38) 

where the superscript "(4)" is to remind us that these matrices are computed in the 
basis of the 4-quasiholes states (14). Consider, for instance, the transformation (35) 

3 what could really change the phase is not the entire electron but the neutral Majorana 
fermion, which is non-local with the quasihole 
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acting on the 4-quasiholes wave functions (14). It changes x — > l/x and 
^(14) (23) so *at 

o(4)x T/ (Q,i) _ (msriu)* ( w , 1 v \ 

*12 T 4qh - / ^(14)(23) ± ^*(13)(24) J ■ 



(13)(24) 



Expressing 77237714 = 77137724* an d taking out y ±y 1 /jc from the denominator we 

(4) 

obtain the exchange matrix in the basis (14) as in Eq. (38). Notice that there is 
no more in this matrix. At this point the Yang-Baxter equations (37) imply 

that if there is a braid-group representation over the 4-quasihole wave functions 
(14) they must be obtained from Eq. (36) by multiplying all exchange matrices 
with e' 51 / 8 . Because is diagonal it could be directly obtained by first fusing 
77i 772 and then interpreting the exchange as 7712 — > e i7r 77i2, which gives the same 
result. Indeed, it follows from Eq. (29) that when r}[ = 772 — 7712, 772 = 772, which 

(4) 

is equivalent to R\J the two functions ^0 an d *Fi in Eq. (29) are multiplied by 
1 and i respectively (note that the expression in the (•■ ■) in Eq. (29) is actually 
independent of 771). 

(4) (4) 

Next, instead of computing R 23 directly it is more convenient to compute R\ 3 , 
following Ref. [38], and then use the identity 

R$=R$R$(r$Y 1 . (39) 



To this end we first apply the coordinate transformation 

n'\ = ?73, 77 2 = 772, 773 = T]i, 77^ = 774, such that 77(3 = e ijr 77i 3 , 



x -> x = 1 -x and l J / ( i 3 )(24) -> VJ/ (13)(24), while ¥(14X23) -> *(12)(34)- Then, using 
the Nayak-Wilczek identity [38] (note the misprints in Eq. (3.8) there) 



(1 -*)¥( 12 )(34) = V i'(13)(24) ~ *¥(14)(23) , as w eU as 



1 + v / 1-x±V1-v / 1-jc=±v / 2V 1 ±V^ 5 (40) 



(4) (4) 

it is not difficult to derive R\$, hence obtain R 23 ' by Eq. (39), i.e., 



R 
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The sign ambiguity in front of the square root in the right-hand side of Eq. (40) 
comes from taking a square root and is not directly linked to the sign under the 
square root on the right. In order to fix this sign we use Eq. (39) and the fact that 
the NOT gate L^Q T Sarma for |x| < 1, that we obtained in Eq. (32) directly in terms 
of the 4-quasiholes wave functions monodromies in Sect. 5, is actually the square 

(4) 

of the exchange matrix , i.e., 

^(4)^(0,1)^(1,0) ( R W) 2 _(0 l\ 

Remark 2 The relevance of the Pfaffian qubit for quantum computation can be 
emphasized once again. Because the monodromy matrix (^23 ) coincides with 

(4) 

the NOT gate in the 4-quasiholes states basis, the matrix R 23 should be identified 
with VNOT, which cannot be implemented in classical information theory [48]. As 
we shall see later this operation is crucial for the construction of the Hadamard 
gate, which is one of the most important single-qubit quantum gate. 

Remark 3 The derivation of the elementary exchange matrices (38)) generating 
the entire two-dimensional representation of the braid group SSa,, over the 4-quasiholes 
Pfaffian wave functions, follows the lines ofRef. [38] where the first row of the ex- 

(4) 

change matrix R\^ has been explicitly computed. These matrices can be obtained 
from the general representation of the braid group £$A,for the Pfaffian state, as de- 
rived in Ref. [39], using the quantum group structure of the parafermion Hall 
states, see Eqs. (138) and (140) there. Nevertheless, the direct and self-contained 
derivation of the exchange matrices, given above, from the analytic properties of 
the Pfaffian wave functions has certain advantages and provides an independent 
check of the results. 

The monodromy transformations, corresponding to the elementary exchanges, i.e., 

(4) 

the squares of R a J a+1 representing a complete counter-clockwise cycle of the parti- 
cle with label a + 1 around that with label a, 

form a subgroup of the braid group called the monodromy group. The monodromy 
group representation is thus generated by the Pauli matrices 0\ and 03, or, al- 
ternatively, by the two elements^ S = i02, R = 0\, such that S 4 — R 2 = I and 
R l SR = S . Therefore the image of the monodromy group is isomorphic to the 
finite non-Abelian group S>^ [49,50], known as the symmetry group of the square, 



4 note that the second Pauli matrix here appears in this context naturally multiplied by i 
and this is in accord with the standard quantum computation conventions [24] 
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which has 8 elements typically given in the two-dimensional representation as 



% = {±I 2 , ±cti, ±io- 2 , ±a 3 }, 



where a ; are the three Pauli matrices. The monodromy group M lx is in general a 
normal subgroup of the braid group SB n and the factor-group 



is isomorphic to the permutation group S^ n . Therefore, the fact that the monodromy 
group representation ^4 = ^4 for the 4-pt functions is finite implies that the 
braid group representation in this case is a finite group whose order could be 
shown to be (Image (^4) | = 96. The order of the representation of the braid group 
^6 is (Image \ = 46080 = 2 6 6!, while that of the monodromy subgroup is 
(Image (^e) I = 32. Similarly, the order of the representation of the braid group 
£$% is |Image(^8) I = 5160960 = 2 7 8!, and that of its monodromy subgroup is 
(Image (^#s) I = 128. These numbers have been obtained by direct enumeration of 
the distinct matrices, produced by multiplying the elementary braid matrices, in 
the corresponding braid-group representation using the Dimino's algorithm [51]. 
In general the image of the representation of the braid group £%2m for n > 3, over 
the 2n-point functions of the Ising model is [52] 



While this may look very nice as a mathematical fact it is fairly disappointing from 
the perspective of topological quantum computation. The point is that our inten- 
tion in TQC is to implement quantum gates by simply exchanging quasiparticles 
positions and the finite braid-group representation implies that we could generate 
only finite number of gates with the Pfaffian qubit. Therefore it cannot be used for 
universal TQC where we would like to efficiently implement any unitary matrix 
(with a given precision). However, the set of quantum gates that could be realized 
by braiding of Pfaffian quasiparticles is known to be a Clifford group, which plays 
a central role in quantum error correction codes, so their topologically protected 
construction is fairly important. Moreover, in the next Section we will try to cir- 
cumvent this restriction by finding a unitary transformation, which does not belong 
to the braid group, that could be used to construct a universal set of gates, and is 
topologically protected. 



Image (^2n) 




for n = even 



for n = odd 
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6.2 The read-out of the single-qubit state in terms of the exchange matrices 

Instead of interpreting the read-out as "taking 773 around 774", as we did in Sects. 3.1 
and 3.2, we could use the exchange matrices (38) to directly "take r/ 3 around r\\ 
and 772". Again the measurement result depends on whether the absolute value of 
the crossratio (1 8) is smaller or bigger than 1 . In the case |x| < 1 , which is the usual 
situation when, e.g., |tJi| > \r\2\ > 1 773 1 > 1 774 1 , as shown on Fig. 7, the operator 
corresponding to the transportation of 773 around 771 and 772 is 




Fig. 7. The read-out for |x| < 1 in terms of the elementary exchange matrices. 

and the result is the same as Eq. (21). If on the other hand \x\ > 1, which can 
be considered as, e.g., (772! < (773!, because using the the CFT symmetry we can 
express x = 773/772. In this case the read-out operation corresponds to Fig. 8, is 




Fig. 8. The read-out for |x| > 1 in terms of the elementary exchange matrices. 



^read-out 12*23 ~ -\ ) \\ Q ) ~ \-\ Oj 

and this is (up to phase) identical with Eq. (24). 



7 Single-qubit gates constructed from elementary exchange matrices 

In this Section we shall implement by quasihole braiding the following single- 
qubit gates: the Hadamard gate H, Pauli gates X, Y, Z, and the phase gate S [24]). 
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In Sect. 8.4 we shall also construct entirely in terms of quasihole braidings the 
Controlled-NOT gate. While not sufficient for universal QC, these gates are known 
to form a Clifford group [40] and could in principle be used for universal quantum 
computation provided that we can create the so-called magic states [40]. The only 
single-qubit gate that we cannot construct directly by braiding is the tt/8 gate T 
[24]). However, instead of the k/S gate one we could use the three-qubit Toffoli 
gate [24]). 

The Hadamard gate H is of central importance for any QC scheme. It is worth 
stressing that the Hadamard gate is the only gate which must be added to a univer- 
sal classical computer (based on the Toffoli gate) in order to make it a universal 
quantum computer [53]. In the TQC with Pfaffian qubits it can be used to create 
special superpositions, called the Bell states (or EPR states) [24], that can be con- 
structed in no other way. In addition H is one of the building blocks of the quantum 
Fourier transform [24]). The Hadamard gate in the Pfaffian TQC scheme can be 
expressed in terms of three elementary braidings of the 4-quasiholes states (14), 
namely [£], 

9 e'i /l A 

H ~ Rj 2 R 13 = #12*23*12 = -i= L J- (41) 




Fig. 9. Braiding diagram for the Hadamard gate (41) and its symbol (on the right) in stan- 
dard quantum-computation notation 

It would be convenient to represent these braidings in braid diagrams, such as 
Fig. (9). The first and the third braids in Fig. (9) are in clockwise direction and 
correspond to the inverse exchanges R^, while the second exchange is in counter- 
clockwise direction and corresponds to R23 . 

The Pauli X gate, known also as the NOT gate, which was first implemented for the 
Pfaffian qubit in Ref. [28], could be executed by two elementary exchanges as 



5 the results in this Section have been originally announced in Ref. [47]; here we give a 
detailed derivation 
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Fig. 10. Braiding diagram for the Pauli X gate and its quantum-computation symbol 

and the corresponding braid diagram is shown on Fig. 10. Similarly, the Y gate, 
which is usually defined in quantum computation literature without the imaginary 
unit, is realized by 4 elementary exchanges as follows 




and the corresponding braid diagram is shown on Fig. 1 1 . The Pauli Z gate [24] can 




Fig. 1 1 . Braiding diagram for the Pauli Y gate 
be realized in two different ways by two elementary exchanges as 

z = r2 \2= R 34 — _] 1 
and the braid diagram for the first of them is shown on Fig. 12. Finally, the phase 




Fig. 12. Braiding diagram for the Pauli Z gate 
gate S, which can also be realized in two different ways by a single elementary 
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exchange, 

S = *12=*34=(J i) 

and the first of them, S = R\2, is shown graphically on Fig. 13. Notice that S 2 = Z 




Fig. 13. Braiding diagram for the phase gate S 

as it should be. 

The only single-qubit gate which cannot be implemented directly in terms of quasi- 
holes braiding is the 7r/8-gate T = diag ^l,e ,jr / 4 J [24] because detT = e ,7r//4 , while 

det yR^a+i) = z f° r a ^ a - I nstea d of T we shall propose to construct the Toffoli 
gate. 



8 Two-qubits construction and two-qubit gates 



In order to realize two qubits, which belong to C 2 , we need at least 4-dimensional 
space. Recalling that the dimension of the excited Pfaffian states with In quasiholes 
at fixed positions [38] is dim j^n = 2" _1 we consider the 6-quasihole states. Before 
we explain how to construct the two-qubit states let us recall that the single qubit 
states can be written as 

|0) = (0+0+0+0+), |1) = (CT + C7_C7 + C7_), 



where we take the convention that the first two a fields determine the state of the 
qubit while the last two guarantee the preservation of the fermion parity, so that 
basically |0) ~ o + o+, while 1 1) ~ o + o , which is in agreement with our definition 
of the qubit because of the fusion rules o+o+ ~ I and o+o ~ \jf. 

The two-qubit basis is defined here by the convention that the first two quasiholes 
form the first qubit while the last two quasiholes form the second qubit 
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m % Va ^ n 6 



Fig. 14. Two qubits, realized by quasiholes with coordinates (t]i,T]2,T]3,T]4) and 
(t/3 , T]4, T]5, Tj6), spanning the needed 4 dimensional space C 2 . The state of qubit 1 is de- 
termined by the quasiholes with positions (t] 1,172), while the state of qubit 2 by the quasi- 
holes with positions (TJs,^) and that is why these two groups were shaded. The state of the 
quasiholes with coordinates (773, T74) depends on both quasihole pairs (t]i, 172) and (175, T]6). 

1 00) = (a + a + c7 + rj + rj + c7 + ), |01) = (rj + a + <7 + c7_rj + rj_) 

1 10) = (o + o o + o o + o + ), |11) = (o+O-O+o+o+O-). (42) 

This is convenient since if we fuse the first two quasiholes this would project to the 
second qubit, while if we fuse the last two quasiholes this would project to the first 
qubit, i.e., 

|aj8) - |/3), |aj8> -> \a). (43) 

Then the third and the fourth quasiholes are fixed by the conservation of the fermion 
parity, i.e., if e,- is the parity of a ep consider the correlation function with yp plugged 
in the middle. This gives 

{<J ei O e2 O ei YFOe 4 O e5 Oe 6 ) = ^1^3 {o ei <J e2 <J e3 O e4 <J e5 <7e 6 ) 

if we move yp to the left or e4e$e( ) (c ei o e2 c e . i c eir O e5 c eb ) if we move it to the right. 
Therefore we obtain the fermion parity rule 

where the last equality follows from ef = 1. Thus we have only 4 independent states 
in the space of 6-quasiholes with fixed positions, which correspond to Eq. (42). 

8.1 Exchange matrices for 6-quasiholes 

The braid group representation over the 6-pt functions is generated by the elemen- 
tary exchanges R^, R%j, R^, R$ and We shall construct them explicitly by 
using the operator-product expansions of the Ising model and the projections to the 
single-qubit states along the lines of Ref. [47]. 
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The construction (42) of the two-qubit states allows to derive the 6-quasiholes ex- 
change matrices in terms of those for the 4 quasiholes. In order to obtain the ex- 
change Rfy we may first fuse 775 — > because the state of the first qubit is inde- 
pendent of 775 and T]6, and use the projections (43). In more detail, we have 



|00) -> {0+0+0+0+), |01) -> (a+a+o+CT_ip) 

110) -> (cr + c7_c7 + c7_), 111) — >■ (<7 + <7_(7 + <7 + y/). (44) 

(4) 

Now the exchange r/i 772 is represented by the action of R\i, which we take 
from Eq. (38), i.e., 



|00) -> (0+0+0+0+) , |01) -> (o+o+o+a_ip) 
110) — >• i (o+<7_o+ct_), 111) — ► i (o+CT-O+o+y) 

11^2 J?1^J?2 



and restoring back the second qubit we obtain in the basis (42) 



R 



(6) 
12 



/l 0\ 

10 

0/0 

\0 // 



i? ( 12 } ®I 2 . 



(45) 



Now let us compute the exchange matrix R^ . Again we can fuse 775 — > r]6 and use 

Eq. (44). The exchange 772 <-> T73 is represented by the action of i? 23 , from Eq. (38), 
i.e., 



|00) -> — = (0+0+0+0+) -i — = (o+0-0+0-), 
m^m V2 V2 



e'4 



e'4 



|01) -> — =(o+o+o+o-\\r)-i — =(o+0-0+o+y), 
m^m V2 v2 



|10) M ^(^-^-)-^(^ + 0- +< 7 + ), 



e ; 4 



e ! 4 



^2^173 V2 V2 



and restoring back again the second qubit according to Eq. (44) we obtain in the 
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basis (42) 



R 



(6) 
23 



e 4 
V2 



( 1 





— / 
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— i 


—i 
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— i 





1 



(4) 
23 



(46) 



Precisely in the same way, by fusing first rji — > r] 2 , we can obtain the 6-quasiholes 
exchange matrices 



? (6) 
c 45 
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(6) 
56 
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h®R$, and 



f (4) 

c 34 ■ 



(47) 



(48) 



The exchange matrix R^ cannot be obtained simply in this way because the quasi- 
holes at T]3 and T74 depend on the states of both qubits so that fusing either (171,172) 
or (r]5,r]e) would change the state of the quasihole pair at (173, 174). As we shall see 
in Sect. 8.2 this entanglement of the two qubits allows us to construct immediately 
the Controlled-Z gate. Notice, however, that must be diagonal and therefore 
could be directly determined by simply using the OPE for the fusion 173 — > 174 alone. 
One way to see this is that if R$ were non-diagonal the exchange of quasiholes at 
rj3 and 774 would create a coherent superposition of the states I|NS) and i//|NS) in 
the Neveu-Schwartz sector (states with even number of o fields acting on the vac- 
uum belong to the NS sector) which would be a violation of the superselection rule 
defined by the chiral fermion parity yp. In contrast, if there are odd number of a's 
to the right of r] a and 77 a+ i acting on the vacuum then the fusion rj a — > T7 a+ i gener- 
ates I|R) and y|R), which are in the Ramond sector where the chiral fermion parity 
is spontaneously broken [33], so that the above states could indeed form coherent 

superpositions. This explains why R^ and R^j- could be non-diagonal, while R^, 

R^4 and R^ have to be diagonal. Thus, using the (neutral part of the) OPE (27), 
we have 



|00) -> 77 34 1/8 (c7 + c7 + c7 + c7 + ), |01) -> r] 3 ^ (o+o+yo+O-) 

T73^I?4 f?3-*»?4 

1 10) -> vi i 3 , A \e + 0-yo + e + ), |11) ^ 77 34 1/8 (a+(7-CT + c7-). (49) 

t?3— *i74 m^ri4 



Therefore, the exchange 773 <-> 174, which simply transforms T734 — > e' n r}i,4, leads to 



34 



|00> 



01) 



10) 



11) 



773^774 



-il -1/8/ \ 

e (0+0+0+0+), 

. — i K 3/8 / v 

?e 'TTJ^ (c7 + C7 + V/(7 + C7_) 

. _;g 3/8/ v 
ze '*Tfc£ (CJ+CJ-V/C7+C7+), 

-i? -I/ 8 / \ 

e ^T7 34 (o+o^o+O-). 



Taking into account that there is another r/^ , coming from the w(l) component of 
the quasihole operator, we find the exchange matrix in the basis (42) to be 



Unlike the other 6 quasiholes exchange matrices, R y 3 ^ is not a factorized tensor 
product of the 4 quasiholes exchange matrices. Instead, there is an additional built- 
in structure in this representation of the braid group, which will eventually allow us 
to construct the Controlled-Z gate. 

It is easy to check that the 6-quasiholes exchange matrices (45), (46), (47), (48) and 
(50) indeed satisfy the Artin relations (37) for the braid group [2], including the 
Yang-Baxter equations. As mentioned before, using the Dimino's algorithm [51] 
we can explicitly obtain the entire group generated by the matrices (45), (46), (50), 
(47) and (48), giving the orders of the representation of the braid group SS^ and its 
monodromy subgroup 

|Image(^ 6 ) I = 46080, |Image(^ 6 ) | = 32. 



8.2 The Controlled-Z gate in terms of 6-quasiholes braidings 

Using the explicit expressions for the 6-pt 7?-matrices, Eqs. (45), (50) and (48), it 
is straight forward to construct the most important two-qubit gates — the CNOT or 
CZ gates in terms of the braid matrices, e.g., 



R 



34 



(I 0\ 
0/00 

0/0 

\0 1/ 



(50) 




/l 0\ 

10 

1 

\0 -I J 



(51) 
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Some insight into the CZ construction may be gained from the identity [24] 

CZ = e'? e'? ZlZ2 e~ ! '? Zl e~'? Z2 , 



where Zi and Z2 are the Z gates acting on qubit 1 and 2 respectively. Because these 
matrices are diagonal and square to I it is not difficult to prove that their exponents 
are actually equal (up to overall phases) to the matrices R^ 2 , respectively, while 
the exponent of Z1Z2 is proportional to the inverse of . 

The braid diagram for the 6-quasihole exchanges corresponding to Eq. (51) is 
shown on Fig. 15. In plotting Fig. 15 we have used that the three i?-matrices en- 
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1 



Fig. 15. The braid diagram for the Controlled-Z gate realized by 3 commuting elemen- 
tary 6-quasiparticle braidings defined in Eq. (51). The symbol on the right is the standard 
quantum-computation notation for CZ. 

tering Eq. (51) are diagonal and therefore commute, which also follows from the 
Artin relations (37), so that the order of the exchanges is not important. Note the 
remarkable simplicity of this realization of the CZ gate — just three elementary ex- 
changes. This is one of the main advantages of the two-qubit construction in terms 
of 6 quasiholes presented in Sect. 8. 



8.3 Single-qubit gates in the two-qubit basis 



Before we continue, it is important to show that we can efficiently express the 
single-qubit gates into the two-qubit basis (42). The point is that the exchange ma- 
trices for 4 quasiholes, which represent the single-qubit operations, belong to the 
braid group ^4 while those for the two-qubit gates are expressed in terms of braid 
matrices from the braid group ^ and the former have different structure from 
the latter. Physically, the embedding of the one-qubit gates into the two-qubit sys- 
tem is non-trivial because the entanglement creates non-local effects between the 
two qubits. Nevertheless, the single-qubit construction in terms of 4 quasiholes ex- 
changes is certainly instructive for the representation of these gates in the two-qubit 
basis. For our purposes it would be convenient to construct these gates explicitly. 
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The Hadamard gate acting on the first qubit can be expressed as 




/l 1 0\ 
10 1 



10-1 
\0 1 0-1/ 



while that acting on the second qubit should be identified with 




1 
1 









1 
1 





1 



1/ 



(52) 



Both Hadamard gates have similar structures to their single-qubit counterparts, yet, 
they are slightly different. This is surprising because adding an additional qubit is 
equivalent to introducing two more strands in the braid diagram and we expect that 
the two straight lines representing a trivial braiding in the new qubit should corre- 
spond to the the unit operator in a tensor product with the single-qubit gate acting 
on the old qubit. The point is that, however, the representation of the braid group 
3§(y realized by the 6-quasihole Pfaffian wave functions naturally appear in a differ- 
ent basis, which is not a factorized tensor product of the representations of SS^ over 
the 4-quasiholes Pfaffian wave functions. This is some kind of topological entan- 
glement which seems to be common for all TQC schemes based on non-Abelian 
anyons realized in FQH systems. 

The phase gates acting on the first and second qubits are respectively 
S 1 =S®I 2 = R { S and S 2 = h®S = B$. 

The embeddings of the other single-qubits into the two-qubit basis follow from 
these of H and S because Z = S 2 and HZH = X. 

8.4 The Controlled-NOT gate 

Now that we know how to construct the Controlled-Z gate, and how to embed the 
single-qubits gates, entirely in terms of 6-quasiholes braidings, the CNOT gate is 
readily computed with the help of the target-qubit Hadamard gate (52), i.e., 



CNOT = H2 CZ H2 = ^56^45^% ^34 R12R45R56 



(\ 0\ 

10 

1 

\0 1 0/ 



(53) 
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We used here that the target-qubit Hadamard gate could be executed by 3 ex- 
changes, according to Eq. (52), as well as the property (R56) 4 = I4 and some of 
the Artin relations (37) for the 6-quasiholes exchange matrices. An equivalent re- 
alization of CNOT, which gives precisely the same result as in Eq. (53), could be 
given in terms of other 7 elementary exchanges 

CNOT = RJR45R34R12R56R45R3I • (54) 




Fig. 16. The braid diagram for the Controlled-NOT gate executed by 7 elementary 6-quasi- 
particle braidings corresponding to Eq. (54). The symbol on the right is the standard quan- 
tum-computation notation for CNOT. 

This is the first known construction of the Controlled-Z and Controlled-NOT gates 
entirely in terms of the braid matrices for 6 quasiholes in the Pfaffian TQC scheme, 
which certainly guarantees the exactness and topological protections of these gates. 

Note that this construction of the CNOT gate is equivalent to the braid realization 
of the Bell matrix of Refs. [54,26]. The algebraic structure behind this Bell matrix, 
when it is used as a universale matrix in the R(T ®T) = (T <g> T)R relations, giving 
rise to an exotic new bialgebra called 503, has been clarified in Ref. [55]. 



8.5 The Bravyi-Kitaev Controlled-Z gate precursor 



For the sake of completeness and comparison we shall also describe the existing 
idea [29] to realize CZ by taking one quasihole, around two other, which suffers 
from the drawback that the resulting transformation has one extra minus sign and 
thus has to be supplemented by external operations in order to produce the CZ gate, 
see below. Consider, e.g., the quasihole with position r\\, from qubit 1 transported 
adiabatically around the two quasiholes, with positions 775 and r/6, of qubit 2 (or, 
equivalently, taking the two quasiholes comprising qubit 2 around one quasihole 
of qubit 1) as shown on Fig. 17. This is obviously equivalent to first taking 77 1 
around 7]^ and then around 775 so that this gate would be just the product of the two 
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Fig. 17. The Controlled-Z gate precursor in terms of the monodromies CZ = Rf 5 Rf 6 . 



corresponding monodromies 



CZ — Rf 5 Ri(, ^ ^12^23 ^34 ^45^56^45^34^23^12 



/l 0\ 

0-100 
10 
\0 0-1/ 



• (55) 



The second equality in Eq. (55) is just an equivalent representation which can be 
readoff from Fig. 17. We give for convenience also the explicit expression for R^ 5 



^15 = ^12 1 ^23 1 ^34 1 ^45^34^23^12 = e' 4 



/ 1\ 

0-10 

10 

\-l 0/ 



which together with Eq. (56) below can be used to compute CZ directly. The braid 
diagram for this realization of the CZ gate precursor is shown on Fig. 18. Now 
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Fig. 18. The braid diagram for the monodromy -based Controlled-Z gate precursor. 

it is obvious that the gate we constructed in Eq. (55) differs from the CZ gate by 
having one more minus sign. The reason is that if the second qubit is in the state 
1 1) the transport of the quasihole at 771 will always produce a minus sign whatever 
the state of the first qubit (note that the two states which are multiplied by —1 
by our gate (55) are exactly |01) and |11)). The idea of Bravyi-Kitaev [29] is to 
split the first qubit into two 1/4-charge states only if this qubit is in the state |1) 
and then move the two quasiholes (at positions 775 and 775 in our case) forming 
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the second qubit around the first qubit in order to execute this gate. This would 
remove the minus sign from the second row of CZ in Eq. (55) and, if successfully 
implemented, should give us a topologically protected Controlled-Z gate . 



8.6 Realization of the Bravyi-Kitaev two-qubit gate #3 



One particular universal set of quantum gates, relevant for TQC with Pfaffian qubits, 
which has been proposed by Bravyi and Kitaev, is [29] 
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The two-qubit gate g2 is identical with our Controlled-Z gate (51) implemented in 
a topologically protected manner by 6-quasiholes braidings. The single-qubit gate 
gi, known also as the tt/8 gate T, has been realized in Ref. [29] as an unprotected 
gate, by bringing together the two quasiholes for a short period of time and then 
pulling them back, in which the exponential topological protection is lost. 

In trying to construct the gate #3 it would be instructive to compute first the mon- 
odromy R^ 6 as taking r\\ around rjg. As is obvious from Fig. 19, this monodromy 
can be expressed in terms of the elementary 6-quasiholes exchanges^] (omitting the 
superscript (6) in the jR-matrices) as 




Fig. 19. The monodromy R 2 l6 executed by taking T]i around % in a counter-clockwise 
direction. 



p2 p— 1 p— 1 p— 1 p— 1 r>2 P P P P 

K 16 — K \2 K 23 K 34 K 45 ^56^45^34^23^12 



/0 1\ 

10 

10 

\1 0/ 



(56) 



note that the clockwise exchanges correspond to the inverse exchange matrices 
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Then the exchange Rk, is just the "square-root" of Eq. (56), i.e., 



^16 = ^i2 1 ^23 1 ^34^45 1 ^56^45^34^23^12- 



(57) 



The direct computation, using the explicit formulas, Eqs. (45), (46), (50), (47) and 
(48) for the elementary R matrices, shows that the two-qubit gate #3 identical with 
R l6 , i.e., 
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The braid diagram for the realization of R^ is shown in Fig. 20. 

^16 = ^'tt"X5^6W3M2 
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Fig. 20. Implementing the two-qubit gate by counter-clockwise exchange of the quasi- 
holes with positions 771 and r]^. 



Remark 4 It is worth stressing that the TQC scheme of Das Sarma et al. [28] is 
essentially based on the monodromy transformations of the multi-quasihole Pfaf- 
fian wave functions. The only exception is the Bravyi-Kitaev gate #3, which was 
proposed to be constructed schematically in Ref. [29] by braiding two among 4 
quasiholes {not by braiding 6 quasiholes as it should be). Note, however, that as a 
two-qubit gate acting on Pfaffian qubits, #3 must be constructed in terms of trans- 
formations of the 6-quasiholes states, which has not been done in Ref. [29]. In 
particular, it is not clear from the 4-quasihole braiding construction of Ref . [29], 
which two quasiholes among the 6 ones must be exchanged in order to obtain the 
gate gi, and whether this is at all possible. The first explicit results proving the use- 
fulness of braidings for topological quantum computation with Pfaffian quasiholes 
have been obtained in Ref. [47] and Eq. (57) is the braiding implementation of ^3. 

Just for reference, and to demonstrate the importance0 of the choice of quasiholes 
to be exchanged, we note that in a similar way the exchange of 772 with 775 gives 



7 note that R25 and R\(, are not equivalent as quantum operations 
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rise to the following matrix 



^25 — ^23 ^34 ^45^34^23 



/ 1 -i\ 

e'"f 1/0 

71 i 1 

\-i 1 / 



5. 7 77ze non-demolition measurement gate 



One of the quantum gates in the universal TQC schemes reviewed in Ref. [29] is 
executed by a non-demolition measurement of the total topological charge of two 
qubits. According to our construction of the two-qubit states this measurement is 
equivalent to the transformation 

|00)->|00), |01> — »• — 101>, 1 10> — ► — 1 10>, 1 1 1> — > |01>. 



Therefore we can identify the non-demolition measurement gate with the mon- 
odromy matrix corresponding to the adiabatic transport of T73 around 774 



K0 : 



/1 o\ 
0-100 
00-10 

\0 1/ 



Z1Z2, 



where Zy = \ R$ J and Z 2 = \ R!g J are the Pauli Z-gates over the first and second 
qubits, respectively. While the non-demolition measurement of the total topological 
charge might happen to be noisy, the implementation of the above quantum gate as 
a 6-quasiholes state monodromy is completely protected by topology. 



8.8 The two-qubit swap gate 



Once we now how to construct the CNOT gate it is straight forward to obtain the 
two-qubit Swap gate in terms of three CNOTs [24] as shown in the top line of 
Fig. 21. Here we shall demonstrate that it is possible to implement the Swap gate 
with 15 elementary exchanges. The bottom line of Fig. 21 shows how to express 
the swapped CNOT in terms of the CZ and Hadamard gates. It turns out that in this 
circuit H is essentially equivalent to so that substituting [H ® I2) — and 
(I2 ® H) ~ 2? 4 5 , as well as using Eq. (51) for CZ, we finally obtain 



Swap ~ (I 2 (g> H) CZ (I 2 ®H){H®li) CZ (H®I 2 ) (h <8> H) CZ (I 2 ® H) 
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- H 



CO 



- H 



H - 



— H 



H - 



H - 



Fig. 21. The two-qubit SWAP gate realized by three CNOT's 

: ^45^34^12^56^45^23^34^12^56^23^45^34^12^56^45 
(\ 0\ 

10 
10 
\0 1/ 



(58) 



9 Universal TQC scheme based on the Hadamard gate H, phase gate 5, 
CNOT and Toffoli gate 



One of the standard universal set of gates, which could be used for universal quan- 
tum computation, includes the Hadamard gate H, the phase gate S, the two-qubit 
CNOT gate and the Toffoli gate, which is a three-qubit Controlled-Controlled-NOT 
(CCNOT) gate [24]. In order to execute three-qubit gates, such as the Toffoli and 
Fredkin gates, we need to consider a system with 8 quasiholes, whose Hilbert sub- 
space of states (for fixed positions of the quasiholes) has dimension 22 _1 = 8. Here 
we shall assume that the third qubit is defined by 2 more quasiholes at positions 777 
and 7]8 as shown on Fig. 22. Then the three-qubit basis can be written in terms of 
the Ising spin fields as follows 



q ubit 1 qubit 2 qubit 3 




»7i ^2 v 3 n 4 v 5 v 6 Vi v s 



Fig. 22. Three qubits constructed from 8 quasiholes 



1 000) = (a + a + a + (7 + c7 + c7 + c7 + c7 + ), |001) = {o + o + o + o o + o + o + o ■_) , 

|010) = (<7 + <7 + <7 + C7_<T + C7_C7 + C7 + ), |011) = ( 0+ <7 + <J+ 0+ 0+ O + O ) , 
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1 100) = (<7 + C7 C7+CT C7 + C7 + C7 + C7 + ), 1 101 > = ( C + O C + (J + C + + + O ) , 

1 1 10) = (<7 + <7 C7+CT+C7+C7 CT + C7 + ), 1 1 1 1) = ( <7 + <7_ <7 + <7_ <7 + <7_ (7+ C7_ ) .(59) 

Now the fermion parity conservation requires that 



<? 3 e 4 = eie 2 e 5 e 6 e 7 es, 



which reduces the number of independent states to from 16 to 8. Although we have 
chosen the quasiholes at T73 and 774 in such a way to preserve the fermion parity, 
any other choice of the positions of quasiholes representing the three qubits would 
be equivalent to Eq. (59) because the braid matrices for the elementary exchanges 
would be just conjugated by an element of the braid group, and the Artin relations 
(37) are invariant under conjugation. 



9.1 Exchange matrices for 8 quasiholes 



Using the fusion rules (27) of the non-Abelian quasiholes we can express the ex- 
change matrices for 8 quasiholes recursively in terms of those for 6 quasiholes as 
follows: 



R { $ = diag(l , 1 , 1 , 1 , i, i, i, i) = R\y ® I 2 , 



(6) 
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= diag(l, l,l,i) 
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(61) 



(62) 



(63) 



(64) 



44 



/ 1 








— i 











o\ 





1 


— i 




















— i 


1 

















—i 








1 


























1 








— i 

















1 


— i 




















— i 


1 
















—i 








1 J 



(65) 



R%> = diag( 1, i, 1, i, 1, i, 1, /) = 1 2 <g> 2? g . (66) 

It is not difficult to check that the exchange matrices (60), (61), (62), (63), (64), 
(65) and (66) satisfy the Artin relations (37) for the braid group SS%. Again, the 
order of the representation of the braid group S$% and its monodromy subgroup can 
be obtained by Dimino's algorithm [51] to be 

|Image(^ 8 )| = 5160960, |Image(^f 8 )| = 128. 



As an illustration of the derivation of the 8-quasiholes exchange matrices, let us 

(8) 

compute the last row of R 67 , i.e., we consider the transformation of the state 1 1 1 1) 
when we exchange T]6 with r]j. Because the state of the second and the third qubits 
is independent of the quasiholes at rji and 772, we could fuse rji — > 772 obtaining in 
this way a 6-quasiholes state whose braiding properties are already known. Indeed, 
using the OPE (27) we find 



1 1 1 1 ) m ^ J ^ ( ¥(12) o+(t?3) a- (t]4)o+ (tj 5 ) o (He) 0+ im) a (ns) ) 



\/ ^J^( <7 -( T ?3)^-(T]4)o- + (T]5)o--(T] 6 )c7 + (T] 7 )c7_(T7 8 ))~|01) 

m^m V 2 y 27723 

where we used the OPE [46] ^(^2)^3(^3) - (277 2 3) _1/2 a- e3 ( T ?3), for 772 -> 773, 
and the identity {o c c + c c+c ) ~ {c + c+c + o c + c ) = |01). It is now obvi- 
ous that the exchange of 776 with 777 in the three-qubit state 1 1 1 1) is equivalent to 

the exchange of the fifth and sixth quasiholes in the state |01) so that, taking k£) 
from Eq. (47), we obtain 



in.) R M 4./^,/ ' H|oo ) + |oi)) 



V2 V 2 V 27723 
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if / 

^= y ^ H( V(i?2) (tj 3 ) <t+ (r] 4 ) o+(i?s) (T?e)^+ M o+(?7 
+ {^(T?2)<7_(773)a + (r74)o- + (T7 5 )o-_(77 6 )o- + (T77)c7_(T78)>) 

- ~7= (-»<CT+(ni)CT_(Tfc)CT + (Tj3)CT_(n4)CT+(T?5)0+(TJ 6 )0+(Tf7)CT + (l?i 

+ (o-+(T?i)^-(^2)a + (7] 3 )c7_(774)(J + (775)a_(77 6 )a + (77 7 )(7_(77 8 ))) 
= ^(-/|100) + |111)), (67) 



(8) 

which exactly reproduces the last row of R 61 ' . In the above derivation we restored 
the 8-quasiholes states using the same OPEs for 772 — ► T73 and 771 — > 772, however in 
reverse, as well as used the identity o (773) cr + (774) ~ c + (773)^(774). 

Remark 5 Due to the specifics of the braid group representation, it may not be 
always possible to represent exactly the single- and two- qubit gates in the three- 
qubit basis (59). Indeed, the 6-quasiholes exchange matrix R^, defined in Eq. (50), 
is not a factorized tensor product of the exchange matrices for 4 quasiholes, rather 
it contains the built-in CZ matrix ( 51 ). Therefore some tensor products, which are 
trivial otherwise, might not be constructed directly in terms of the elementary ex- 
change matrices for 8 quasiholes. One consequence of this peculiarity is that some 
one-qubit and two-qubit gates would be easier realizable in the three-qubit basis 
(59) in terms of elementary 8-quasiholes exchange matrices, however, up to a pair 
of extra minus signs, see Sect. 9.2. While the exact construction of the one-qubit 
and two-qubit gates would require more work, their simplified versions might be 
sufficient in most cases. 



9.2 Embedding of one-qubit and two-qubit gates into a three-qubit system 



The three one-qubit phase gates are directly expressed as single elementary 8- 
quasiholes exchange matrices, i.e., 



The first one-qubit Hadamard gate can be constructed exactly in terms of the ex- 
change matrices for 8 quasiholes by (skipping the superscript "(8)" of R) 
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Hi = H ® I4 ~ Rn R 23 R n 
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, (68) 



while the second Hadamard gate could be constructed as follows 
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The third Hadamard gate could also be reproduced up to some swapping by 



H 3 = I 4 ® H ~ 2? 78 2? 45 2? 56 R 67 R 56 R 45 R 1S 
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(70) 



It might be useful to give also the realization of the single-qubit NOT gates: X\ = 
X <g> I 4 = Rj 3 , X 2 = h ® X <8> l 2 = #45 and X 3 = I 4 ®X = R^Rjj. 

We should stress again that the topological entanglement between the qubits men- 
tioned in Sect. 8.3 leads to serious difficulties for efficient embedding of the one- 
qubit and two-qubit gates in three-qubit systems. For example, the NOT gate X 3 
acting on the third qubit has a different structure than just a tensor product of the 
exchange matrix producing the single-qubit X gate. Similarly, the Hadamard gates 
acting on the first and second qubits have slightly different structures from their 
single-qubit counterpart H, while that acting on the third qubit cannot even be ob- 
tained exactly with the same number of elementary exchanges as H. This seems 
to be a common problem arising in all TQC schemes using non-Abelian anyons 
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in FQH systems, whose general solution is still missing. Moreover, it appears that 
the two-qubit Controlled-NOT gates in a three-qubit systems cannot be directly 
constructed in the three-qubit basis (59) because of the topological entanglement 
between the two qubits and the third one. This requires more work and will be 
reported elsewhere. Just for reference, we give a simple implementation of a three- 
qubit operation which is very close to CNOT2 = I2 <8> CNOT 



CNOT 2 ~ I2 <8> CNOT ~ /?7 2 1 /?5 6 SWAP 2 /?36/?45SWAP 2 
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(71) 



where 7? 36 = R^R^R^R^R^ and the gate SWAP2 is defined below. 

The two-qubit SWAP gates can be simply realized by braiding in the three-qubit 
basis (59) (up to overall phases and pairs of extra minus signs) by 



SWAPj ~ SWAP ® I 2 ~ R 4S R^R^R 4S R 73 R^R^R 2 3R4sRx R 3i ^45 
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(72) 



SWAP 2 ~ I 2 ® SWAP ~ R 61 R^R 5 ^R 61 
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The Swap gates turn out to be very important because they can be used to con- 
struct gates acting on one of the qubits in terms of similar gates acting on another 
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qubit. For example, the CNOT gate acting on the first and the third qubits can be 
expressed in terms of the CNOT acting on the first and the second qubit plus two 
gates SWAP2. The extra minus signs appearing in some of the three-qubit opera- 



-0- 



Fig. 23. The CNOT acting on qubits 1 and 3 expressed in terms of CNOT® I2, which is the 
CNOT gate acting on qubits 1 and 2, and two gates I2 <8> SWAP. 

tions are not an innocent thing, because they may have different properties from 
the standard gates. However, it appears that in many cases these simplified gates, 
which are much simpler to construct, can be used instead of the standard once. 

The three-qubit Toffoli gate [24] can be constructed in terms of the Controlled-S 
gate and CNOT like in Ref. [47] using the relation between the Toffoli gate and 
the Controlled-Controlled-Z gate or by a braid-group based Controlled-Controlled- 
Z gate precursor Ref. [47], which must be supplemented by the Bravyi-Kitaev 
construction [47]. 



10 Discussion 



In this paper we explicitly implemented all single-qubit gates in the Pfaffian TQC 
scheme, except for the n/S one, in terms of 4-quasihole braidings, as well as the 
two-qubit Controlled-Z and CNOT gates in terms of 6-quasihole braidings. These 
gates, which are known to form a Clifford group, are realized in a completely topo- 
logically protected way because of the topological nature of the braid operations 
in the FQH liquids. This work is an extension of the topological quantum com- 
putation scheme of Ref. [28] using pairs of Pfaffian quasiholes localized on anti- 
dots to construct elementary qubits and execute logical NOT on them. While the 
original TQC scheme of Ref. [28] used only monodromy transformations to real- 
ize quantum gates, we, for the first time, exploited explicitly quasihole braiding in 
the Pfaffian FQH state to construct the single-qubit Hadamard gate H, the phase 
gate S and the CNOT gate. Although the Gottesmann-Knill theorem says that any 
circuit based only on the Clifford group gates could be efficiently simulated on 
a (probabilistic) classical computer these gates do play a crucial role in quantum 
computation, especially in the error-correcting algorithms [24] . 

Due to the topological entanglement between the separate qubits realized by non- 
Abelian anyons in FQH systems some difficulties arise when trying to embed the 
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one-qubit and two-qubit gates into systems with more qubits. This makes the em- 
bedding of Clifford gates non-trivial and requires more work. 

For implementing three-qubit gates such as the Toffoli and Fredkin gates [24], in the 
Pfaffian TQC scheme, we considered Pfaffian wave functions with 8-quasiholes, in 
which case the topological degeneracy of the space of correlation functions is 8 
[12,38]). We derived explicitly the braid matrices for the elementary 8-quasiholes 
exchanges, which serve as building blocks for constructing all three-qubit gates. 
More work in this direction, including eventually the topologically protected con- 
struction of the Toffoli gate would be reported elsewhere. 

To conclude, let us make some remarks about the possible observation of the non- 
Abelian statistics. We believe that in order to observe the Pfaffian phase at v = 5/2 
one should perform the experiment at temperature below 15 mK. The point is that 
there might exists another incompressible Abelian phase, which was called the Ex- 
tended Pfaffian (EPf) state [33], that could also be realized at v = 5/2. Perhaps the 
most observable difference between the two phases is in the electric charge of their 
quasiparticles: 1/4 for the Pfaffian and 1 /2 for the EPf. The EPf phase was obtained 
mathematically by a local chiral algebra extension of the Pfaffian state and satisfies 
all conditions necessary for an incompressible quantum Hall state [33]. The moti- 
vation for introducing this new state is that there is a persistent unexplained kink 
around T c = 15 mK observed in the thermal activation experiment [56] showing 
two different slops that presumably correspond to two different gaps below and 
above the critical temperature. Analyzing the edge states CFT it has been proposed 
in Ref. [33] a possible explanation of the kink in terms of a finite temperature 
two-step phase transition between the Pfaffian and the EPf state involving an in- 
termediate compressible state of composite fermions. Here is a brief description 
of the process (see Sect. 9 in Ref. [33] for more details): at low temperature the 
system is definitely in the Pfaffian phase, as the numerical calculations suggest. As 
temperature increases to about 1/2 of the Pfaffian-phase gap, which was estimated 
to be about 33 mK (for the sample of Ref. [56]), the system becomes more and 
more compressible (look at the behavior of the free energy on the edge, Fig. 5 in 
Ref. [33]) leading to a Il-nd order phase transition to the compressible state of com- 
posite fermions (which has the same topological structure like the Pfaffian state but 
having at the same time the Z2 symmetry of the EPf state that is broken sponta- 
neously in the Pfaffian phase). Immediately after that, as temperature continues to 
increase, there is a I-st order phase transition to the EPF state which is expected to 
have a higher gap than the Pfaffian state. 
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A Binomial series expansion of the 4-pt function and analytic continuation 

In this appendix we shall give some details about the analytic continuation of the 
function (19), which might be necessary for the understanding of the results of 
Sect. 3. Using the standard complex analysis notion [57] of a branching point as 
a multi-valued isolated singular point, we consider a punctured neighborhood of 
the branching point, denoted as U' in which we would like to continue the element 
([/,/) of the function (19) from the simply-connected sub-domain U C U' along 
any path. For example, for the branching point at z = we can choose 

U' = {z\0< \ z \< 1}, U = {z \ \z-l/2\ < 1/2}, 
as shown on Fig. A. 1 . Then the analytic continuation along the contour 




Fig. A.l. (Color online). Domains and contours for the different values of \z\ used for the 
analytic continuation and binomial series expansion 
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7o 



z=\**, 0<t<2n 



changes the sign in front of the inner root but not the one of the outer root because 
7o does not encircle z = 1 . This can be verified directly by using the (fractional 
powers) Laurent expansion, which in this case can be obtained by the binomial 
series expansion 

(!+,)« =f r ^ + 1 ) J M<1 



applied for a = 1/2. Using the defining T function property T(z) = (z — l)T(z— 1), 
we get 



1^1 1 



That is why yj\±yjz — > a/1 Ty/z, for |z| < 1, when continuing z — > e 2?n z along 
any contour in C/', which is homo topic to (A.l). 

For the branching point at z = 1 , on the other hand, we consider the domains V' and 
V shown again on Fig. A.l defined by 

V' = {z\0<\z-l\ < 1}, V = {z \ \z-l/2\ < 1/2}. 



Then the analytic continuation of the element (V,f) of the function (19) from V to 
V' along any contour in V' homotopic to 



7i 



z=l + -e u , 0<t<2n 
2 



does not change the inner root sign because the point z = is outside the contour. It 
only changes the sign of the outer root whenever the inner root sign is "— ". Indeed, 

representing y/z = a/1 + (z-1)~ (l + ^ - + ^-), for |z- 1| < 1, and 

using again the binomial expansion in V' we get 



z-\ 



1 + ^=^1 l + ^--—(z-l) 2 + ---) (A.3) 



1 -*=^\}--r+m ( *- 1 r + -"h zeV ' (A - 4) 
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The appearance of \/z — \ in the right-hand side of Eq. (A.4) but not in Eq. (A. 3) 
implies that when z encircles the point 1 the function \/l — \fz acquires one minus 



sign, while the function J\ + ^fz is single- valued. 

Recall that for a contour passing through the points z = or z = 1 it is not possible 
to make analytic continuation. 

Finally for \z\ > 1 we use W' = {z | 1 < \z\ < °°}, W = {z \ \z — 2\ < 1}, and the 
continuation along the contour 
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jz z = ie\ 0<t<2n} 



changes both the sing of the inner root and that of the outer root when the inner 
sign is "— ". Because in this case 1 1 / y^z\ < 1 the binomial expansion with respect 
to\/y/z gives 



1±Vz= ^i 1± _i__^ ± _i_... 1 , 1<lz ,< 



and that explicitly shows that z = °° is a branching point of order 4. 

That is how we conclude that the general contour y, which i shown on Fig. A.l, 
corresponding to the read-out transformation (17) is homotopic to Jq or fo U 7i 
depending on the value of |x|, i.e., 



7o f° r |*| < 1 

7o U 71 for \x\ > 1 



which explains once again the analytic continuation results obtained in Sect. 3. 
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